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Abstract 

A systematic density functional theory (DFT) +U study is conducted to investigate the electron 
correlation and spin-orbit coupling (SOC) effects in US3 and USe3. Our calculations reveal that 
inclusion of the U term is essential to get energy band gaps for them, indicating the strong cor- 
relation effects for uranium 5/ electrons. Taking consideration of the SOC effect results in small 
reduction on the electronic band gaps of US3 and USe3, but largely changes the energy band shapes 
around the Fermi energy. As a result, US3 has a direct band gap while USe3 has an indirect one. 
Our calculations predict that both US3 and USe3 are antiferromagnetic insulators, in agreement 
with corresponding experimental results. Based on our DFT+t7 calculations, we systematically 
present the ground-state electronic, mechanical and Raman properties for US3 and USe3. 
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I. INTRODUCTION 



Actinide based materials possess interesting physical behaviors due to the existence of 
strongly-correlated 5/ electrons and have attracted extensive attentions 1-6 . Different from 
actinide oxides which have been widely studied to reveal the detailed configurations and 
correlation effects of the 5/ electrons^— , none of the actinide chalcogenides has ever re- 
ceived comparable concerns. Here we take US3 and USe3 as two representatives to study 
the electron correlation and spin-orbit coupling (SOC) effects of 5/ electrons in actinide 
chalcogenides. The other reason for us to investigate their electronic structures is that they 
employ the layered MX 3 structure (with M to be a metal element, and X to be S, Se, or 
Te), which can be used as models for studying electronic behaviors in 1-dimensional (ID) 
systems. The MX3 structure belongs to the space group of P2i/m, with its atomic organi- 
zations depicted in Fig. 1(a). In the monoclinic lattice, the top two chalcogen (S , Se or 
Te) atoms form a tightly bound pair suggesting that the MX3 compounds may be regarded 
as MX(X 2 ) 3 -. The top X-X pair together with an underneath M and a more underneath X 
atom form a triangular unit, which repeats along the b lattice direction forming a ID chain, 
and two such triangle units form a unit cell of the MX 3 compounds in the (a,c) plane. In this 
way, the MX3 compounds are always considered as ID materials. For example, the metallic 
ZrTe3 and NbSe3 exhibiting charge-density-wave transitions have been thoroughly investi- 
gated in angle-resolved photoelectron spectroscopy (ARPES) and optical experiments^ - —. 
In the present paper, by contrast, we address the electronic structures of US3 and USe3, two 
semiconducting members of the MX 3 family. 

In studies of actinide based materials, one has to be careful for the electron correlation and 
SOC effects of the actinide 5/ electrons. As an example, for actinide dioxides conventional 
density functional theory (DFT) schemes that apply the local density approximation (LDA) 
or the generalized gradient approximation (GGA) underestimate the strong on-site Coulomb 
repulsion of the 5/ electrons and consequently fail to capture the insulating properties' 1 ^. 
Several approaches, the LDA/GGA+£7, the hybrid density functional of (Heyd, Scuseria, 
and Enzerhof) HSE, the self-interaction corrected local spin-density (SIC-LSD), and the 
Dynamical Mean-Field Theory (DMFT), have been developed to correct the pure LDA/GGA 
failures in calculations of actinide materials. Among them the effective modification of pure 
DFT by LDA/GGA+J7 formalisms has been widely used in theoretical studies of UO^^ 
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and PuO^r^ - — >2S. The obtained structural parameters as well as the electronic structure 
and phonon dispersion curves^2& accord well with experiments. In our present work, the 
GGA+[/ schemes due to Dudarev et a/.— - — are employed to study the electron correlation 
and spin-orbit coupling effects in US3 and USe3, as well as the two materials' electronic, 
mechanical, and Raman properties. 

The rest of the paper is organized as follows. In Sec. II the computational method is 
briefly described. In Sec. Ill we present the results of the physical properties of US3 and 
USe3, and discuss the electron correlation and SOC effects of the uranium 5/ electrons. 
Finally in Sec. IV, we close our paper with a summary of our main results. 

II. CALCULATION METHOD 

Our total-energy calculations are carried out by employing the plane-wave basis pseu- 
dopotential method as implemented in Vienna ab initio simulation package (VASP}2£. The 
exchange and correlation effects are described with the GGA approximation in the Perdew- 
Burke-Ernzerhof (PBE) forrn^^. The projected augmented wave (PAW) method of Blochl^ 
is employed with the frozen-core approximation. Electron wave function is expanded in plane 
waves up to a cutoff energy of 400 eV and all atoms are fully relaxed until the Hellmann- 
Feynman forces on them are less than 0.01 eV/A. A 7x9x1 Monkhorst-Pack— k-point 
mesh is employed for integration over the Brillouin zone of US3 and USe 3 . The uranium 
6s 2 6p 6 5f 3 6d l 7s 2 , sulphur 3s 2 3p 6 , and selenium 4s 2 4p 6 electrons are treated as valence elec- 
trons. Noncollinear calculations are used when considering the spin-orbital coupling effects. 
The strong on-site Coulomb repulsions among the localized uranium 5/ electrons are de- 
scribed by using the formalism formulated by Dudarev et al.— — . In this scheme, the total 
GGA energy functional is of the form 

E GGA+ u = E GGA + $^[Trp CT - Tr(pV)], (1) 

cr 

where p u is the density matrix of / states with spin a, while U and J are the spherically 
averaged screened Coulomb energy and the exchange energy, respectively. 

In this paper, the Coulomb U is treated as a variable, while the exchange energy is set 
to be a constant J=0.51 eV. This value of J is in the ball park of the commonly accepted 
one for uranium compounds 2 ^ and close to the theoretically predicted value of 0.54 eV in 
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UO2 42 . Since only the difference between U and J is significant, we will henceforth label 
them as one single parameter U e g=U-J, while keeping in mind that the nonzero J has been 
used during calculations. 

The phonon frequencies at the Gamma point for US3 and USe3 are calculated by using 
the density functional perturbation theory (DFPT). And their Raman-active frequencies 
are obtained through symmetry analysis on the corresponding vibration modes. Before 
DFPT calculations, the lattice constants and atomic positions of US3 and USe3 are further 
optimized using a denser fc-point mesh of 9x13x3, and a finer force convergence criteria of 
0.001 eV/A. A 2x2x1 supercell is subsequently used for DFPT calculations. Since both 
US3 and USe3 are stacked along the c direction through van der Waals interactions, and 
with relatively large lattice constants, we do not extend the supercell along the c direction. 
During DFPT calculations on the 2x2x1 supercell, the /c-point mesh is set to be 5x7x5, 
and the energy convergence criteria is set to be le-6 eV. 

III. RESULTS AND DISCUSSION 

Experimentally, the magnetic orderings of uranium trichalcogenides were studied ever 
since 1961 3 i>^~— . The relatively large atomic distances between neighboring U atoms in 
US 3 and USe 3 suggest that magnetic orderings may occur in these compounds through 
super-exchange interactions of uranium 5/ electrons via the S or Se ions 3 -^. And different 
measurements indicate that both US3 and USe3 crystals undergo AFM transitions at very 
low temperatures, with the magnetic moments of the two uranium atoms within each unit 
cell opposite to each other— In the present GGA+U study, we consider the nonmagnetic 
(NM), ferromagnetic (FM), and antiferromagnetic (AFM) phases for each choice of the U e g 
value and then determine the lowest-energy state by a subsequent total-energy comparison 
of these three phases. Our calculated electronic energies for different magnetic phases of US3 
and USe3 are all listed in Table I. One can see that within GGA formalism or the GGA+U 
formalism with a too small U e s value, the FM state is more stable for both US3 and USe3, 
in contradiction with experimental results. Therefore, the electron correlation effect has to 
be accounted for to correctly describe ground-state US3 and USe3. In the discussions that 
follow, we will confine our reports to the AFM solutions for US3 and USe3. 

The experimentally determined lattice parameters of US3 and USe3 are (a=5. 37 A, £>=3.96 
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A, c=9.94 A) and (a=5.65 A, 6=4.06 A, c=10.47 A) respectively^^. Our calculated 
lattice constants with different U e s values are summarized in Table II, together with the 
experimental results. We can see that different from what we found for actinide dioxides^ 1 ^, 
the lattice constants of actinide trichalcogenides do not change monotonically with the value 
of the U e g parameter. Besides, the lattice constants obtained in GGA calculations are 
obviously too small in comparison with corresponding experimental results, for both US3 
and USe 3 . For US 3 , the value of U e g has to be as large as 6 eV to get reasonable lattice 
constants compared with experimental results. Differently, a U g r value of 4 eV is enough to 
get reasonable lattice constants for USe 3 , and changing U e s from 4 to 6 eV has negligible 
effects on the obtained lattice constants. 

Besides of the prominent changes in the atomic structure parameters, the most dramatic 
improvement brought by the GGA+U formalism when compared to the GGA results is in 
the description of electronic structure properties. For this, we have investigated the band 
structures in AFM phases for US3 and USe3 aiming at revealing the fundamental influences 
of considering the on-site Coulomb interaction. Figures 2(a) and 2(b) show the obtained 
local density of states (LDOS) for the SI, S3, and U atoms in US 3 , and Sel, Se3, and U atoms 
in USe3 respectively. Without accounting for the on-site Coulomb repulsion, one can see 
that the GGA calculations predict an incorrect metallic ground state by nonzero occupation 
of uranium electrons at the Fermi energy (Ef). When switching on the U e s parameter, as 
shown both in Figs. 2(a) and 2(b), the uranium electronic states in US3 and USe3 begin 
to split at E{ and tend to form two peaks with the gap of AE. The amplitude of this gap 
increases with increasing U e g. Previous electric resistivity measurements pointed out that 
USe3 did not conduct at room temperatures^, thus a band gap at Ef should be contained 
in its electronic structure. We can see from Fig. 2(b) that an energy band gap appears 
only when Ugg is larger than or equal to 4 eV for USe3. For US3, the value of U e Q has to be 
enlarged to 6 eV to open the band gap at the Fermi energy. These results also imply that the 
electron correlation strength might be different in US3 and USe3. From the last figure in Figs. 
2(a) and 2(b), we can see that the sulphur (selenium) and uranium electronic states overlap 
with each other and contribute equally to the valence band maximum (VBM) of US3 (USes), 
while the conduction band minimum (CBM) is composed of uranium electronic states. The 
orbital mixing between sulphur (selenium) and uranium electronic states below the Fermi 
energy indicates that there are covalent interactions between the sulphur (selenium) and 
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uranium atoms in US3 (USes). 

To further analyze the orbital-resolved electronic structures, we calculate the projected 
density of states (PDOS) for the X-np (n=3 (X=S) or 4 (X=Se)), U-5/ and XJ-6p electronic 
states in UX 3 . Figures 3(a)-3(d) show the obtained PDOS for US3 in different calculations, 
while the corresponding PDOS results for USe3 are shown in Figs. 3(e)-3(h). The value 
of U e g is chosen to be 6 eV in all GGA+U and GGA+£/+SOC calculations. One can see 
that without considering the electron correlation effect, US3 and USe3 are both wrongly 
predicted as metallic materials, which contradicts with experimental results. After adding 
the U parameter to describe the strong on-site energy of U-5/ electrons, the energy gaps 
are opened for US3 and USe3. We can see from Figs. 3(d) and 3(h) that further considering 
SOC has little influences on the electronic states around the Fermi energy. For both US3 
and USe3, the SOC effect lies in the deep energy level, causing an energy splitting of the 
U-6p states. From the PDOS results, we can also see clear difference between the XI- and 
X3-np (n=3 (X=S) or 4(X=Se)) electronic states. The PDOS for the electronic states of the 
X2 atom is very similar to that of the XI atom and thus is not shown here. The different 
electronic state distribution of the XI (X2) and X3 atoms proves the theory of considering 
UX 3 as UX(X 2 ), and indicates the strong covalent bondings between the XI and X2 atoms. 

To analyze more carefully the SOC effects in US3 and USe3, we recalculate the electronic 
structures of them with higher resolutions, with the reduced Gaussian smearing width of 
only 0.02 eV and a PDOS resolution of 0.01 eV. Figures 4(a) and 4(b) show the obtained 
PDOS results for uranium 5/ electrons of US3 and USe3, by using the GGA method, while 
Figs. 5(a) and 5(b) show the obtained energy band structures of US3 and USe3 by using 
the GGA+U method, respectively. All the GGA+SOC and GGA+[/+SOC calculations are 
noncollinear, with both SOC and orbital polarizations considered. And the results before 
and after considering the SOC effects are shown in solid and dotted lines respectively. From 
the PDOS results obtained by using the GGA method, we can see that pure inclusion of SOC 
to the GGA method does not lead to energy band gaps for US3 and USe3. This discovery is 
different from the situation of CoO, where SOC can solely open a band gap 51 . 

From the band structures of US3 shown in Fig. 5(a), we can see that both the lowest un- 
occupied conduction and highest occupied valence bands distribute along the T-Z direction. 
For the lowest conduction band, the SOC effect causes an energy downshift of 0.05 eV along 
the T-Z direction. Contrarily for its highest valence band, the SOC effect causes an energy 
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upshift of 0.03 eV at the T point and an upshift of 0.07 eV at the Z point. As a result, the 
energy band gap of US3 changes from indirect (T — ;-Z) into direct type (Z— >-Z) after con- 
sidering the SOC effect, and the gap value reduces from 1.59 to be 1.48 eV. Different from 
US3, the lowest conduction and highest valence bands of USe3 distributes along the K^-B 
and B-r-Y directions respectively. For the lowest conduction band along the K^-B direction, 
the SOC effect causes an energy downshift of 0.04 eV at the Ki point and a downshift of 
0.02 eV at the B point. Correspondingly the CBM of USe3 moves from B to the Ki point. 
However, since the VBM of USe3 is at some point along the B-r-Y direction instead of at 
any high-symmetry k points, the SOC effect does not change the indirect band gap type. 
Overall the band gap of USe3 changes from 1.30 to be 1.23 eV after considering the SOC 
effect. Our studies reveal that the energy band gap of USe3 is smaller than that of US3. 

After systematically presenting the electronic structure results for US3 and USe3, we 
now turn to their mechanical properties. The elastic constants are obtained by solving the 
eigenvalues of their Hessian matrix, which is calculated based on the Hooke's law and small 
position changes on independent S, Se and U atoms. The obtained elastic constants of US3 
and USe3 in different calculations by using the GGA and GGA+f/ methods are both listed 
in Table III. We can see that the elastic constants calculated by using the GGA formalism for 
US3 and USe3 do not satisfy the stability criteria" for monoclinic structures that the Cn, 
C22, C33, £44? C55, C 66 should be all positive. This result further proves that conventional 
GGA formalism fails for describing the uranium trichalcogenides. Contrarily, the calculated 
elastic constants by using the GGA+U method for both US3 and USe3 satisfy the above 
stability criteria, as well as the other six criteria for monoclinic structures", proving the 
stable existence of US 3 and USe 3 , and their strong electron correlation effects. 

Based on the elastic constants of US3 and USe3, we further calculate the Voigt and Reuss 
bounds on their bulk (By, Br) and shear moduli (Gy, Gr)'. And the bulk and shear 
moduli of US3 and USe3 can be estimated by B=(By+B^) /2, and G=(Gy+Gn)/2. And 
the Poisson's ratios can be calculated by u=(3B-2G)/(QB+2G). The obtained mechanical 
properties for US3 and USe 3 are all listed in Table IV. We can see that the Voigt and Beuss 
bounds obviously differ from each other for bulk and shear moduli. This result reflects the 
structural anisotropy of US 3 and USe 3 . Besides, the bulk and shear moduli are both found 
to be larger for US3 than for USe3. The Poisson's ratio is calculated to be 0.19 and 0.22 for 
US3 and USe 3 respectively. 
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Until now, there are no experimental reports on the mechanical properties of US3 or 
USe3. Hence our theoretical values can be used as references for further investigations or 
industrial applications of them. At another side, the Raman properties of the MX 3 group 
materials have been systematically studied for a long time, and some of the characteristic 
peaks of US3 and USe 3 were successfully observed in experiments^. Therefore, here we 
further calculate the Raman properties of US3 and USe3. According to group theory and 
the symmetry of them, both US3 and USe3 are predicted to have 12 Raman-active vibrational 
modes^i^: 8 A g and 4 B g modes respectively. The symbols A g and B g represent for two 
different vibration symmetries respectively: vibrations inside and out of the mirror plane. 
Specifically, A g corresponds to the vibrations where dy is always while B g corresponds to 
the vibrations where dy is not 0. 

Our Raman results of US3 and USe3 are obtained by symmetry analysis based on the 
vibrational modes calculated within the GGA+U method, where SOC effects are consid- 
ered. Figures 6(a) and 6(b) show our determined Raman frequencies for US3 and USe3 
respectively, together with the experimental results by G. Nouvel et al.—. One can see that 
by using the GGA formalism for both US3 and USe3, the Raman frequencies are totally 
confused compared with the experimental results. At another side, the GGA+U results on 
the Raman frequencies are clearly reasonable in two obvious aspects: i) there exist a high- 
est single Raman frequency for both US3 and USe3, corresponding to the molecular-like X 2 
mode vibrations; ii) the highest Raman frequency of US3 is much larger than that of USe3, 
indicating that the S-S interaction is stronger than Se-Se. When comparing all the Raman 
frequencies, there seems to be a frequency shift of about 25~35 cm -1 between our GGA+U 
and the experimental results. 

IV. CONCLUSION 

By using the GGA and GGA+U methods, we have systematically studied the electronic, 
mechanical, and Raman properties of US3 and USe3, aiming to reveal the underlying elec- 
tron correlation and SOC effects. By comparing the calculated lattice constants of US3 and 
USe3 with corresponding experimental results, and monitoring their electronic structures 
with different U e s parameters, we conclude that a U c s value of 6 eV is needed to get rea- 
sonable lattice constants, and the right insulating properties for both US3 and USe3. After 
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considering SOC effects, the deep U-6p band will be split into two separate energy bands 
in both US 3 and USe 3 . The SOC effects also changes the energy band shapes around the 
Fermi energies and lead to the fact that US 3 is a direct band gap while USe 3 is an indirect 
band gap insulator. Because of the SOC effects, the energy band gaps of US 3 and USe 3 
reduce by 0.11 and 0.07 eV respectively. Based on the obtained antiferromagnetic ground 
states of US 3 and USe 3 , we then systematically give out their electrical, mechanical, and 
Raman properties. The VBM and CBM of US 3 (USe 3 ) are found to be contributed by 3p-5f 
(4p-5/) hybrid and 5/ (5/) electronic states respectively. The elastic constants of US 3 and 
USe 3 are found to satisfy the stability criteria, and their Poisson's ratios are calculated to 
be 0.19 and 0.22 respectively. For the Raman properties, we theoretically repeat the two 
important experimental observations that the S-S or Se-Se dimer vibrations have the largest 
frequencies in US 3 and USe 3 , and the frequency of the S-S dimer vibration is larger than 
that of the Se-Se dimer vibration. 
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TABLE I: Electronic energies per unit cell of US3 and USe3 in the nonmagnetic (NM), ferromagnetic 
(FM), and antioferromagnetic (AFM) states by using different calculation methods. All values are 
in units of eV. 
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-50.85 


-50.85 


-44.50 


-47.26 


-47.35 


GGA+U {U eS = 


=5 


eV) 


-46.55 


-50.74 


-50.74 


-43.95 


-46.79 


-47.26 


GGA+U {U eS = 


=6 


eV) 


-48.39 


-51.68 


-51.69 


-43.48 


-46.52 


-46.85 
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TABLE II: Lattice constants of US3 and USe-3 obtained by using different methods. All values are 
in units of A. 



"Reference 
^Reference 



47 
48 





us 3 


USe 3 


Methods 








a b c 


a b c 



GGA 




5.11 3.83 9.25 5.50 3.99 10.22 


GGA+C/ (U cS ~- 


= 1 


eV) 5.12 3.82 9.26 5.52 3.95 10.30 


GGA+C/ (U cS = 


=2 


eV) 5.13 3.82 9.28 5.52 3.94 10.42 


GGA+C/ (U eS = 


=3 


eV) 5.09 3.74 9.61 5.50 3.87 10.98 


GGA+C/ (U cS = 


=4 


eV) 5.09 3.75 9.65 5.71 4.02 10.90 


GGA+C/ (U cS ~- 


=5 


eV) 5.06 3.82 9.61 5.73 4.04 10.52 


GGA+C/ (U cS ~- 


=6 


eV) 5.44 3.92 9.76 5.73 4.04 10.60 


Exp. a ' fc 




5.37 3.96 9.94 5.65 4.06 10.47 
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TABLE III: Elastic constants of US 3 and USe 3 obtained in GGA and GGA+f7+S0C calculations. 
The U e fi value is chosen to be 6 eV in both GGA+£7+S0C calculations. All data are in units of 
GPa. 

Cn C22 C33 C44 C55 Cqq C12 C13 C23 C45 Ci6 C26 C36 

US 3 (GGA) -654 -80 -4132 452 -332 -273 -923 -2953 -1766 32 -511 -640 -522 

US 3 (GGA+[/+SOC) 963 1055 437 313 299 169 196 132 181 -6 -102 -55 -74 
USe 3 (GGA) 946 583 -107 363 309 7 -37 -349 -65 3 171 24 198 

USe 3 (GGA+[/+SOC) 951 938 496 257 244 267 195 194 216 -13 -19 -20 -33 
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TABLE IV: Mechanical properties of US3 and USe3 including the bulk and shear moduli, their 
Voigt and Reuss bounds, and the Poisson's ratio calculated by using the GGA+U method. The 
U e g value is chosen to be 6 eV. All moduli data are in units of GPa. 



G\ Gr B G v 
US 3 386 299 286 245 342 266 0.19 
USe 3 400 365 272 257 382 265 0.22 
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List of captions 



Fig.l (Color online), (a) The atomic structures of MX 3 , where blue and dark green balls 
representing for metal (M) and chalcogen (X) atoms respectively. The monoclinic lattices 
are depicted by the dashed lines, (b) Depiction of the Brillouin Zone for the monoclinic 
MX 3 lattice. 



Fig. 2 (Color online). The local density of states (LDOS) for the SI, S3, and U atoms in 
US 3 (a) and Sel, Se2, and U atoms in USe 3 (b) by using the GGA and GGA+£7 methods 
with U e g ranging from 1 to 6 eV. The SI (Sel) and S3 (Se3) atoms correspond to XI and 
X3 atoms depicted in Fig. 1(a). The Fermi energies are denoted by dashed lines. 



Fig. 3 (Color online), (a)-(d) The projected density of states (PDOS) for the 3p electronic 
states of SI and S3 atoms, and 5/, 6p electronic states of the U atom in US 3 . (e)-(h) The 
PDOS for the 4p electronic states of Sel and Se3 atoms, and 5/, 6p electronic states of 
the U atom in USe 3 . The SI (Sel) and S3 (Se3) atoms correspond to XI and X3 atoms 
depicted in Fig. 1(a). The Fermi energies are denoted by dashed lines. The values of U c h 
are chosen to be 6 eV in all GGA+U and GGA+?7+SOC calculations. 



Fig. 4 (Color online). The projected density of states for the uranium 5/ electronic states 
in US 3 (a) and USe 3 (b), by using the GGA method. The Fermi energies are set to be 
zero and denoted by dashed lines. The results before and after considering the spin-orbit 
coupling effects are shown in red solid and blue dotted lines respectively. 

Fig. 5 (Color online). The electronic energy band structures for US 3 (a) and USe 3 
(b) by using the GGA+U method. The values of t/ e ff are chosen to be 6 eV. The Fermi 
energy and different fc-points are denoted by the dashed lines. The results before and af- 
ter considering the spin-orbit coupling effects are shown in solid and dotted lines respectively. 



Fig. 6 (Color online). The calculated frequencies at the T point for Raman-active 
vibrational modes of US 3 (a) and USe 3 (b) by using the GGA and GGA+?7 methods, 



together with the experimental results listed in Ref. |3lJ. The values of U e ^ are chosen to be 
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6 eV. The symmetry for each vibration mode are also presented. 
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